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Abstract 



'sj" . We examine time reversal invariance in the /^-decays of B(2 



and '^Li(2+) to ^Be (2+) in detail, with particular attention to final 
state interactions of the two a particles from the decay of ^Be, and 
of the electron (positron) with the daughter nucleus. An R-matrix 

^ \ formulation is used, with the initial state described by a shell model. 

^ ' The R-matrix parameters are obtained by fitting the a — a scattering 

phase shifts and allowed /3-decay rates. Because the nuclear final 
^ ■ state interaction effects on the time reversal test come from second 

^ . forbidden /3-decays, they are small and can be minimized by a suitable 

choice of kinematic conditions. The e^-nucleus Coulomb interaction 
induced T-violation effects depend on the first forbidden operators. 



1 Introduction 

Tests of time reversal invariance (TRI) remain very important because it is 
only in the K^-K^ system that CP (or T) invariance has been found not to 

*Present address: Physics Department, Fudan University, Shanghai 200433, P. R. 
China. 



hold to about 2 x 10"^ 0. Until TRl violations (TRIV) have been found 
in other systems, it is difficult to distinguish between various theoretical 
proposals [2-7] that have been put forward to explain the experimental results 
obtained in the decays of neutral kaons. 

Nuclear tests of TRI that have been carried out include detailed balance 
[i] {ciT < 5 X 10~^), reciprocity [^ (aj- < 10^^), polarized hyperon decay [|10|, 



TT| , nuclear 7-decay [|T2| and nuclear /3-decay [jl3| (ay < 10~^) experiments. 
The most accurate test of TRI, combined with parity non-conservation is the 
lack of neutron and atomic electric dipole moments |1|, |I^ . For a review of 
experiments and theory, we refer to Ref . |Tl], |T5[ . 

Tests of TRI in /?-decay are made difficult by the necessity of measuring a 
triple correlation or the polarization of the emitted electron. The triple cor- 
relation {J) ■pe^Vu requires a polarized parent state and the measurement of 
the recoil nucleus. The mass eight system is amenable to a somewhat simpler 
test, because the daughter of the decays from ^Li or ^B, namely ^Be*, is un- 
stable and breaks up into two charged particles ^He) which can be detected 
more easily than the recoil. However, a question arises immediately: Does 
the final state interaction amongst the final e^ and the two a particles emit- 
ted in the decay of ^Be* spoil the TRI test? In this work we investigate this 
question in detail. We find that true tests of TRI which depend on measure- 
ments of one (or both) recoil a particles are only spoiled by competing second 
order forbidden /3-decays and by the e^-nucleus Coulomb interactions. By 
a proper choice of kinematics, even this small spoiling can be minimized. 
The primary decay is from a 2"^ state of ^Li or ^B to the lowest 2"'" state of 
^Be*. "Pseudo-TRI violation" (PTRIV) for the strong interaction effects. 



which mimics TRIV through the introduction of phases [0 occur primarily 
through the interference of the main 2+ — > 2"'" matrix elements with those 
of the 2+ -> 0+ and 2+ ^ 4+ (J = 0+ and 4+) of ^Be; the PTRIV effects 
for the e^-nucleus Coulomb interaction depend on the Coulomb phase shifts 
and the 2+ -^ 2+ transition operators. 

A complete ffist-principle theoretical study is immensely complex. We 
treat the decay as a two-step process (see Fig. 1) 



^Li(^B) ^ ^Be* + e-(e^ 



v[v) 



-^ a + a + e (e^) + z/(z/), 

and use an R-matrix formalism in which the initial state is described by a 
shell model wave function and the ^Be* resonant states are saturated by shell 
model states. The parameters of the R-matrix are constrained by physical 
considerations and are adjusted to reproduce a-a scattering phase shifts 
in the relevant energy region. With these parameters, it is found that the 
unpolarized /3-decay spectrum is reproduced with two additional parameters. 
The TRIV and PTRIV correlation observables are studied in this framework. 
In the next section, we examine the TRI tests in the mass 8 system. In 
Section 3, we develop the R-matrix approximation to be used. In Section 
4, we study the PTRIV due to hadronic interactions; in Section 5 we study 
that due to the e^-nucleus Coulomb interactions, and finally in Section 6 we 
summarize our results. 

2 Time reversal invariance tests in the p— 
decay processes of the A=8 System 

The questions related to tests of TRI in the weak and electromagnetic pro- 
cesses have been studied by many authors. We restrict ourselves here to 
the mass 8 system, for which an experimental study has been proposed [p!7|] . 
The initial 2"*" ^Li or *B ground states decay by /3-emission to the (resonant) 
states of *Be* with J^ = 0"*", 2"*", 4+, . . ., which break up into two a particles. 
The graphical representation, together with definitions of some kinematic 
variables are given in Fig. 1. The final 2"*" resonant states of ^Be* are our 
main interest. Those states with J^ = 0^ and 4+ can interfere with the 2"*" 
state to produce PTRIV signals due to final (2a) state interactions. These 
states and interactions should be taken into account in the study of TRIV 
effects. The resonant states of ^Be* are treated in section 3. In the present 
section, emphasis is given to presenting the analytic results and the shell 
model computation of the hadronic weak transition amplitudes. 



a) The decay rates for the unpolarized case 



The differential decay rate for tfie /5-decay of an unpolarized nucleus can be 
expressed as 



dW 



2n5{E, -Ef~te- t,)Tr 



HiY&Hw 



d^, 



(2.1) 



where dC, is an infinitesimal phase space element of the final state, 6 is the 
final state projection operator, H^ is the weak Hamiltonian responsible for 
the weak charged current nuclear reactions and the trace is over the initial 
and final state (spin + linear momentum) subspaces. We write 

1 



Tr 



Hy/QHw 



-Glcos%W''''L^,. 



(2.2) 



When the leptonic polarizations are not detected, W^'^ and L'^'^ are 

W"" = Tr [r\~ci)eJ''i-q)] , (2.3) 

L^"" = Tr [rt(q)P(q)] , (2.4) 

where J^(q) is the Fourier transform of the hadronic weak current operator 
J'^(x, t = 0), /^(q) is Fourier transform of the matrix element of the leptonic 
weak current operators with q = —kg — k,y.Q Eq. \i.l\ can be written as 

dW = 2nS{Ei-Ef-ee-eu) 
^Glcos% 



Ws^L,, + l^i+)^L(+)^ + W^i-)^L(;)^ + W^L^] d^ 



(2.5) 

where i,j = 1. . .3 and irreducible tensor components of the leptonic ten- 
sors Lss, . . . , Lyy are given in the following subsection. Wss and W^^ contain 
contributions of the time (labeled by "s" ) and space (labeled by "f " ) compo- 
nents of the hadronic weak current operators. W^^' contains contributions of 
the interference terms between the time and space components of the same 
current operator, where "+" ("-") denotes the symmetric (anti-symmetric) 
combination with respect to the time and space indices of the corresponding 
operator. 

^■^The Coulomb interactions between the charged leptons and nuclei that are ignored 
here will be investigated in Section 5. 



b) Kinematic functions and hadronic response func- 
tions 

In the decay of *B and ^Li, the daughter ^Be* breaks up into two a particles. 
The resonant state of the two a particles is of the following form 

|-;k),^ = 47r E iZ ^'^e-''■^fY]^.u,mhm^ (2-6) 

J^=0+,2+,4+ Mf=-Jf 

where "-" denotes an incoming boundary condition, k is a unit vector in the 
direction of the relative momentum of the two a particles, Sjj^ is the strong 
and electromagnetic phase shift of the Jf partial wave and Yjm is a spherical 
harmonic. 



Using the trace formula given in \TE], various terms in Eq. ^.5| can be 
expressed in terms of a set of irreducible hadronic response functions, namely, 

WssL^s = {47rfY.^''''~''e'^'-''''"f^Y.Cy'Rss{cr)KUcr), (2.7) 
JfJ'f 

(2.8) 






JfJ'f "P^ 



(2.9) 



where 



CifJ'f 



:J'f\\y^\Jf) 



Jr 



a 



with X = y/2x + 1. The kinematic functions are 



KJa] 



A-K 



Ki^\ap) = 47r[n(k)®[i;(q)®LW], 
K,{apr) = 47r[n(k)®[F^(q)®Lg, 



(2.10) 



(2.11) 
(2.12) 
(2.13) 



where [0/^ ® ^/aljm denotes the Clebsch-Gordan couphng of (pij^mi and 4>i2m2 
to total angular momentum j and magnetic quantum number m. The total 
number of different kinematic functions are 34; 22 of them are T-even and 
12 of them are T-odd. 

It can be shown that the leading terms of the hadronic response functions 
R^stK^^p) and Rvv{(ypT) are of order 0[(i?fi;)^] (it is 0[{RtiY] for i?,,(a)), 
where i? is a typical radius of the parent/daughter nucleus; thus at low 
momentum transfer, terms with large p can be neglected since Rk « 1. 
If one keeps terms up to order 0[k^R'^], the possible kinematic functions 
with 0" < 2 together with their properties under parity and time reversal 
transformations are listed in Table |I|. The T-odd correlation observables 
interested in this paper are P-odd. There are four of them which are given 
in the lower right box of Table |l]. 

The leptonic tensors in case of not observing the polarization of the 
charged leptons are found to be 



Lss 

lL 



[e^ey + ke ■ k,,] , 



8 [eeK + ^u^e] , 

Q 

-j= [SCeC^ - ke ■ k,,] 

±8v/2[e,ke-eek,], 
2. 



l^pl^fy ~r l^iyi^p 



-Kp • K,, 



(2.14) 
(2.15) 
(2.16) 

(2.17) 

(2.18) 
(2.19) 



The hadronic dynamic response functions Rss, • • • , Rm can be expressed 
in terms of the reduced matrix elements of the multipole operators [|T9| 
Cj, . . . , Mj of the hadronic charged weak currents. The analytic form of 
these relations can be found in Appendix A. The T-odd response functions, 
which are of interest in this work, have the following explicit expressions. 
Contribution from the transition to the J^ = 2"*" states: 



RiVm) 



,2 

'5^ 



i^ImCoMl* + i-KJ-ReC^El* 



'h 



(2.20) 



i?(;)(21) = -tyhmCoE;-t'^^hmCoLl-t^V7ImLlC',* 

-i^JllmElC^* + ilJllmLoC;, (2.21) 




i?„„(221) = i^\l^ImL\El* -i^ImL^El+i-^j—ImLlEl* 

i-\—ImL^E^*+i—ReL^,M^* (2.22) 

5V35 ^ ^ 10 ^ ^ ' ^ ' 



i?,,(212) = -ty^ImLlM*-t^^^ImE',M*+i^^hmL,Mt 

+t^^ReE!E; + ^l-^ReE!^ - t-l^RellE;. (2.23) 

Here, symbols with a superscript "5" indicate that they are contributions of 
the hadronic axial vector current operator; others are contributions of the 
hadronic vector current operator. Cj or Cj are reduced matrix elements of 
the multipole operators of the weak charge density operators. Lj . . . Mj are 
matrix elements of the space components of the weak current operators. Lj or 
L J are reduced matrix elements of the multipole operators of the longitudinal 
part of the weak current operators. Ej or Ej are reduced matrix elements of 
the electric part of the transverse components of the weak current multipole 
operators and Mj or Mj are the reduced matrix elements of the magnetic 
part of the transverse components of the same operators. 
Contribution from the interference terms between J^ = 2"*" and J^ = 0"*" states: 

Rit\22) = ^CoMt-^^E!c; + ^C,E;, (2.24) 

R^A2i) = -lf-CoE; + l^^LoC;-lf-c,L;, (2.25) 

R..{221) = -^LoE;-^LlMt, (2.26) 



10V5 " ^ 5V 5 



i?,,(212) = . LoMl*--\-ElEl 




+^LIe;, (2.27) 

where those reduced matrix elements that are complex conjugated correspond 
to transitions to the final O"*" state and all the rest correspond to the transition 
to the final 2"*" state. 
Contribution from the interference terms between J^ = 2"*" and J^ = 4"*" states: 

i?i:'(22) = icoAff + iylsfC^- - ^CiB;, (2.28) 

/2 1 

i?,,(221) = -y^LoE; + —LlMt (2.30) 

lU 15 

fi„„(212) = l^L„Mr + A^B.=B- + i^BfL,- 

-^LlEI, (2.31) 

where those reduced matrix elements that are complex conjugated correspond 
to transitions to the final 4+ state and all the rest correspond to the transition 
to the final 2"'" state. 

The T-odd observables corresponding to these interference response func- 
tions depend on their imaginary parts. 

The T-even response functions can also be written in the same expanded 
form as those given above for the T-odd ones. They will not be presented 
here; their expressions can be deduced by using the compact expressions 
given in Appendix A. 

c) The differential decay date in terms of the weak 
transition ampUtudes 

The kinematic functions can be written in Cartesian form. The differential 
decay rate Eq. ^3| in the 2+ resonant region is 

dW = dW^2^ + dW^2^ + dW^o^ + dW^l> + . . . , (2.32) 



9 



where dW22 is a contribution to the differential decay rate even under time 
reversal and dWjj, are contributions to differential decay rates odd under 
time reversal. The subindices denote the contributing angular momentum 
(J) of the partial waves of the final two a particles. The relatively small 
terms that originate from the interference between the 2"'" and the 0"*" and 
4"*" final states are included for the T-odd observables only since they are 
important contributions to the PTRIV that originate from the final state 
interaction. The contributions from pure 0^ and 4"*" states to the TRIV 
signal are negligible in the energy region where the lowest 2"*" state resonates; 
they are not included in the differential decay rate. 

When TUf; is neglected, the T-even differential decay rate has the following 
generic form 



dwi$ 



dWo 



Mn4%-k. + 4°M(ke-k.)' 



+T2(k) : [{nf^kSe + R?kk + Rf%k 

+ (i?f keke + Rf^keK + R^^KK) k ' K]} , 

where the leading decay rate dWo is 



(2.33) 



dWr 



16 

= {2nf5iE,- Ef - e,- e,)—eee,Glcos%F{Z,e,)\gA\^\A2\''d^, 

5 



(2.34) 



and the non-relativistic one body Gamow-Teller matrix element A2 is defined 

as 



A, 



{2^m 



20 



0" 



2+) 



(2.35) 



with yjim = [Yi (g) e] 



Jm^ 



1 and F{Z, e^) the Coulomb function for the 



charged leptons. In Eq. p.34| , the tensorial contraction between the second 
rank tensor T2(k) and a pair of vectors AB is defined as 



T2(k) : AB 
T2{k\j 



T2{k\,AiB 



J' 



niKj Q^ij- 



(2.36) 
(2.37) 
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The neglect of rrie is justifiable since the energy released in the transitions 



3-10 MeV) is much larger than me. 
Three T-odd differential rates are 



(0) 



dWi^ 

dwS^ 



dWoT2{k) : [4'^keke X k, + i?f k.k^ X k 
dWoT^ik) : [i?f iceke X k, + i?i?k,ke X k 
dWoT^ik) : 



Rukeke X k^ + R^i^KK X K 



where we also have neglected rrie- 

The Cartesian response functions R\ 



, 3) and R, 



(2) 



(2.38) 
(2.39) 
(2.40) 

-,•••,12) 

• • , -^vv • 

They are given in Appendix B. At low momentum transfers, we expand these 
momentum dependent multipole operators of the hadronic charged weak cur- 
rents in terms of their static multipoles, up to order of 0{k,/M) with M the 
mass of a nucleon, and/or 0{R'^k'^) with i? ~ 1 — 2 f m the size of a typical 
nucleus. As a first step, we take the non-relativistic one body approximation 
for the hadronic charged weak current operators by retaining only the leading 
term in an expansion. The justification of such a step can be understood 



are related to those in the spherical harmonics basis |T3] , namely R, 



M 



following the discussions of Ref . [20] . The results are well known |T9[ . They 
are reproduced in Appendix C for completeness. 

The result of the expansion of the hadronic charged weak currents opera- 
tor in the non-relativistic one body approximation and the definition of the 
related static multipole operators are given in Tables || and |. The symbol 
( ||. . .|| ) represents the reduced matrix elements. 

It is useful to introduce the following quantities 



^* = ^r 



A2 



with Aj defined in Table |3|, and 



/i 



/. 



M 




(2.41) 

(2.42) 
(2.43) 
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9 a 



(2.44) 
(2.45) 



where the isospin indices "(±)" for the single nucleon form factors are sup- 
pressed. 

The Cartesian form of the T-odd response functions -R7 , . . . , R\2 ^-re 



R'i 



(2) 



Timfl + Imfu - Imfx 



R. 



(2) 



R, 



(2) 



R 



(2) 
10 



2M 



+ 2M^ 1 -W T^'"^^'^ ^ ~'"'^ ^ 15 VT'""^ ^ r5v/7 



Imr]io 
(2.46) 



A-e. 



2M 
A2 / 1 / 7 



±/m/,^ + ImfM + /m/T 



MM 14 V 15'"^^^^^ ^ ll'^^^ ^ m y'"^^ ^ IE7^'"^^° 
Aee (blY 3v^ ^1/10, +4 



M2 






feA 
"M2 



sm(52 - 5o) W^i?e/ir76(2+,0+) ± ^i?er79(2+, 0+) j 



ImrjiQ 

(2.47) 



+ ^^^^(^2 - <5o) I ^y|i?e/ir^6(2+,0+) ± ^i?er^9(2+, 0+) ) , (2.48) 
-^s^n((52 - (5o) [ ^W^i?e/i^6(2+,0+) ± ^Rer),{2^ ,0' 



12 






^2.49) 



d(2) 



'-^^^^2 - 6,) rJ^RefM2\^^) ± f^J^Rev,{2\A^ 



■^szn{6, - 6,) I ± J^i?eA%(2+,4+) ± 1 J^i^er/9(2+,4^ 



i? 



A 



2 3 



2 2 



12 = ^5m(52-(54) -W-i?e/ir^6(2+,4+)±-W-i?er/9(2+,4 



where 



Aee 



szn{S2 - 6,) I ^^^i?e/ir/6(2+,4+) ± ^^'^Rer],(2^,4^ 



+ '^,sin{52 - 6,) I Aj|i?e/ir^6(2+,4+) ± lJlRe^,{2\A' 



:2.50) 



(2.51) 



^^j+||^{±V2y^||2+) 



^6(2+, J^ 



M 



2 i=l 



5:(2+||rW3^io-a||2+) 
X:(J+||rWr23^22-a||2+) 



^9(2+, J+) 



M' 



2i=l 



E(2'"lk^^^3^io-^||2+) 



1=1 



(2.52) 



(2.53) 



(2.54) 



with J^ = O"*", 4"*". Terms in Eqs. p.48| - p751| that are proportional to cos{. . .) 
have been dropped, since they are much smaller than those terms propor- 
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tional to sin{. . .) due to the fact they are multiphed by T-odd hadronic 
response functions which are expected to be much smaller than the corre- 
sponding T-even components even if the Hamiltonian of the system is not 
invariant under time reversal. The Fermi matrix element Ai is assumed to 
be zero. 

d) Shell model calculation of the matrix elements of 
the static multipole operators in the A=8 system 

The intermediate states of the /5-decay processes are ^Be* resonances. Here, 
we shall consider a shell model computation of the reduced matrix elements 
of various static multipole operators that are used in the next section where 
the decay of these shell model states are discussed and where a comparison 
with experiments is made. 

We shall restrict ourselves to a Ip shell model space. Using the effec- 
tive interaction developed by Cohen and Kurath [^(2BME 8-16), the shell 
model Hamiltonian can be diagonalized [^. The Ip configuration mixing 
shell model calculation with the Cohen-Kurath potential (2BME 8-16) is 
used to calculate the J^ = 2+ levels of *Be* (see Fig. 2); these are stable 
shell model states that lie at E=3.55(0), 13.37(0), 16.19(1), 17.45(0),.. .(the 
numbers in parentheses are the isospin for the corresponding states and the 
unit is in MeV). These theoretically determined energy levels lie reasonably 
close to the experimental resonant state energies given in Fig. 3. The shell 
model state at E=3.55 corresponds to the E=3.04 state and that at E=16.19 
corresponds to the E=16.63 state in Fig. 2. The shell model state at E=13.37 
possesses similar one particle parentage coefficient in its wave function ||2^ as 



the shell model state located at E=16.92, so although it is relatively far from 
the state at E=16.92, it should be identified with that state. The E=17.45 
state corresponds to some higher excitation state in the J^ = 2+ state series; 
it is not considered in this paper. 

The doubly reduced matrix element of an one body transition operator 
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has the following form W^, E4 



<J,;T^\0^\J2;T2> = Y.^j,T{a'a) < a'\\0%\\a>, (2.55) 

aa' 

where a, a' enumerate single nucleon shell model states. In the Ip shell model 
subspace, a = j = (1/2,3/2), with j the single nucleon total angular mo- 
mentum. The density matrix '?/;j.r(a', a) is defined in Ref. [^; it depends on 
the shell model Hamiltonian in the Ip shell subspace. The reduced matrix 
elements of various static multipole operators defined in Table § are 

A = ±-^Tr^^A„ (2.56) 

where "T" denotes transpose and tp is the corresponding density matrix. 

By using harmonic oscillator radial wave functions [^, the one nucleon 
matrix elements corresponding to those defined in Table ^ are 

(2.58) 
A2 = ^^^( ^^ f%], (2.59) 



< 


It > 


^ 


^4n 


< 


It > 




^An 


< 


It > 


s/At: 


< 


It > 


1 


^An 


< 


It > 


1 


^An 


< 


It > 



-1 -V5, 



-2 V2 ' 

-V2 -VTo. 



(2.60) 
(2.61) 
(2.62) 
(2.63) 
(2.64) 
As = ^4^( A °Vm6)2, (2.65) 

(2.66) 



A, 



Ac 



A 



10 



T / 5 /2 

<It> -T 



V47r 


< It > 


\/An 


< It > 


y/An 


< It > 



Att 



10 



10 cr /5 



10 

5 /l" 
2 V 3 



"■V 2 



5 2 
2 V 3 

3 



(M6)^ 



{Mbf, 



^V 2 



(M6)^ 








-^Nte 



{Mh)\ 
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(2.68) 
(2.69) 
(2.70) 
(2.71) 
(2.72) 
(2.73) 
(2.74) 
(2.75) 



where the matrices Aj (i=l,. . . ,10) act on the j = 1/2 and j = 3/2 spaces 
with the 11 (first) element of the matrices corresponding to the diagonal 
matrix element for j = 1/2. The single nucleon isospin factor < It >= a/3- 
b is the oscillator parameter for the shell model orbits. 



3 Semi-phenomenological R— matrix treatment 
for the A = 8 system 

a) R— Matrix theory for /3— decay processes 

Following Appendix D, the T-matrix for the /9-decay processes is found to 
be 



T, 



fi 



E 



Rn{E) is complex, namely 

Rn{E) 



^2a,n(^2a)(^;e^/|ij"^o|0^) 
E2a — En — Rn{E2a) 



Dn{E)~ih{E), 



(3.1) 



(3.2) 
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where the real functions Dn{E) and In{E) are the level shift and the half 
widthp^ of the nth level, respectively. r2Q-,Ti(-E') is the vertex function that 
connects the nth level resonance to the two a states. The result given in 
Eq. |3.1| is formally exact, provided | n) is a complete set in the Hilbert space 
of interest. All effects of AV are represented by the functions Rn{E) and 
Tn{E). We shall not use a complete set of (n |; we only need a few of the most 
important shell model states, as shown in Appendix D. In addition, only the 
(strong and electromagnetic) interaction effects at short distances (r < 5fm) 
are included in Rn{E) and r„(£'). The residual effects are written as 

Tf,= Y: {■■■) + STf,. (3-3) 

shell model 

where residual part STji will be specified in the following. 

The imaginary part of Rn{E) is related to the magnitude of Tn{E) through 
the optical theorem p6|]. In the particular energy range (0 < E < 16 MeV) 



considered here, the only open channel into which the shell model states | n) 
can decay is the 2a one. The optical theorem requires that 



IniE) = 7^p(E)|^2„,„(^)|^ (3.4) 



piE) = ^M^^fM^, (3.5) 

where we have assumed the center of mass (cm.) frame for ^Be* and used a 
non-relativistic approximation 

with k the relative momentum between the two a particles. 

At low energies, the R-matrix parameters In{E) and r2a,n{E) are inde- 
pendent of the details of the interaction potential AV. Their asymptotic 
behavior can be derived from Eqs. p. 7] , and p.8| . Namely, if A^ has a short 
range (the long range Coulomb interaction between two a particles is treated 

■^■^The width of the nth level is defined as the value of In{E) at the energy where 
E ^ En - Dn{E) ^ 0. 
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separately in this study) and is less singular tliarQ 1/r^ "^ (e << 1) as r — i> 0, 
then 

|2«) '"'J^" (a|k|)^, (3.7) 

I ipn) ~ constant, (3.8) 

where a is some scale factor of dimension length that is roughly the size of 
/W and J is the total angular momentum of the state with two asymptotic 
a particles as r ^ oo. By combining Eqs. |^, ^^ and p.7| , one gets 



T2aAE) ""-^ Ei, (3.9) 

where use has been made of the non-relativistic spectrum for the a particles 
Eq. U. 

From Eq. |3.4|, we get 



In{E) ""-^ E-^+i (3.10) 

At low energies, Dn{E) and In{E) are expected to be smooth functions of E. 
They can be expanded in a power series of E with the expansion coefficients 
treated as parameters. These parameters are determined from experimental 
data for a-a scattering phase shifts and allowed /5-decay rates in the A = 8 
system. 

We shall require that D'^{E) is analytic in the energy region of interest. 
The denominator E — E^ — D:^{E) can therefore be written as 

E-Ei- Di{E) = E- E^ - a{E - E^'-^f - b{E - E^'-^f + . .(3,11) 

with E^'"^ the experimental resonant energy of the nth level of total angular 
momentum J. Terms of higher power in E — E^'"^ represented by ". . . " are 



■^■^Albeit the definition of r can be ambiguous if the a particle is composite and the 
overlap between the two a particles is strong, the singularity of the effective AV can be 
smoother than the one extracted from AV at larger distance where the overlap between 
these particles are negligible due to the compositeness of the particles. Therefore the 
statement based on point particle pairs given here is also expected to be true for composite 
particles. 
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expected to be unimportant around resonant energies. The parameters a and 
h are to be determined in the following. 

Next, let's consider the vertex function. From the low energy behavior 
given by Eq. ^^ ^n{.E) can be written as 

Vi{E) = NoEi{l-cE + ...), (3.12) 

with ". . ." representing higher order terms in powers of E that can be ignored 

at low energies, and A'"o, c are coefficients. The half width function li{E) is 

related to T'^{E) through the optical theorem, Eq. |3.4| . In the following, we 

shall express the magnitude of VJ^{E), in terms of Ini^) ^^^ parameterize 
rJ 



KXE) as 




IniE) = /rj-^o h^O 77 ^0777' (3-13) 



where 1°'"^, and c are coefficients to be determined. The vertex function 
T^{E) is assumed to be real function of E on the real ii^-axis. 

b) a— a scattering 

The R-matrix theory can also be applied to the scattering of two a particles. 
The collision of two a particles can be viewed as proceeding through various 
resonant levels of the A = 8 system, plus potential scattering between the two 
a particles at larger distances where the overlap between the two a particles 
is small. The rigorous relationship between the R-matrix parameters in the 
/?-decay processes and the resonant scattering cross section of the a particles 
is indirect. Thus, a more explicit ansatz has to be made in order to relate 
the /3-decay and the a-a scattering processes in a way that is convenient for 
the phenomenological applications. 

We assume that the scattering phase shifts of the a particles due to in- 
teractions at short distances (the resonant scattering), 6j (E), and at large 
distances (the potential scattering), 5j (E), are additive. The resonant 
part of the T-matrix for the collision of the two a-particles in a state of 
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total angular moinentum J can be related to phase shifts generated by the 
resonant states and by the potential scattering as 

tj{E) = --A^e^^^(^)sm<5,(i?), (3.14) 

5j{E) = 5^r'\E) + 5f"'\E), (3.15) 

We shall assume further that the vertex function in the R-matrix does not 
contain any contribution from the potential scattering at large distances, and 
that the potential scattering is provided solely by the long range Coulomb 
interaction between the two a particles. The phase shift due to the potential 
scattering is expressed as 

tan6f°'\E) = -Fj{rk)/Gj{rk), (3.16) 

where Fj{x) and Gj{x) are the regular and irregular Coulomb functions 
[p^ for the two a-particle system, and r is the distance between the two 
a particles where their relative radial wave functions vanish when AV = 0. 
It follows that there are infinite number of discrete values of r that lead to 



the same potential phase shift given by Eq. 3.16 ; all of them are energy 
dependent. The energy dependence of r is expected to be smaller when it is 
close to the origin, especially when the interaction between the two a particles 
is stiff enough at small distances that it can be treated as a hard core. The 
energy dependence of r is thus expected to be minimized by choosing it fairly 
close to the origin P^ . 

The T-matrix corresponding to the interactions that are responsible for 
the resonant states (see Fig. 5) in the A = 8 system can be expressed in 

■^■■^Since we can extrapolate the Coulomb wave function given by a combination of Fj 
and Gj at large distances all the way to points close to the origin where the strong 
interaction dominates, the value of r is not necessarily restricted to the region where the 
Coulomb interaction is most important. In fact, any value of r that gives the same energy 
dependence for the phase shifts is acceptable. The meaning of r as a measure of the range 
of the potential is meaningful only if the potential indeed contains a stiff piece (e.g., hard 
core) . In this case, a value of r that is located near the hard core should be almost energy 
independent. Given the energy dependence of Sj , it must be emphasized that, in the 
region where the strong interaction is important, the zeros of the true wave function of 
the system are not located at the ones satisfying Eq. 3.16| 
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terms of the R-matrix parameters, namely 

_ V^ \92a,n\'^)\ /o 1 7\ 

^ E - E;i - DJ{E) + tw^^iEY ^ ' ^ 

where it is assumed that the level shift function D:^{E) is the same as that in 
the R-matrix. The relationship between g^{E) and w:^{E) can be obtained 
by unitarity considerations of the a-a T-matrix; on the other hand, the 
relationship between w:^{E) and li{E) is not obvious. r;[(i?) is the half 
width of the corresponding resonant level or the inverse life time of that 
level. The decay of the level to two a particles can be pictured as that of 
a prepared state, namely as the eigenstate of the shell model Hamiltonian 
Hsheih tunneling through the potential "barrier" A^. From the time reversal 
invariance consideration, it follows that in case of resonant a-a scattering, 
the time needed for the a particles to tunnel into the potential "barrier" 
AV to form a well prepared resonant state (corresponding to a shell model 
eigenstate) is the same as the time it takes the same state to tunnel through 
the potential "barrier" AV on its way out. Thus the time delay due to 
the a particle staying in a resonant state in a-a scattering is twice as long 
as the time for the corresponding shell model eigenstate to decay near a 
resonant energy of that state. This qualitative argument is formally, albeit 



not explicitly, derived in Ref. |26|. From these considerations, it is assumed 
that 

wi{E) = \li{E). (3.18) 

As will be shown in the following, this relation leads to an excellent synthesis 
between the /?-decay a spectra and a-a scattering within the R-matrix 
formalism and without the need of any "intruder" states. 

Since a two spin zero particle state with definite total angular momentum 
J has only one phase shift at a given energy, the phase for each of the terms 
in the above sum has to be the same, which allows us to write Eq. |3.17| as 



tan5^r\E) = -V- '^^^^\,,^, . (3.19) 

' ^ ^ ^E-Ei-Di{E) ^ > 
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Together with Eq. p.l6| , this equation yields 



Before going on to further specifications of li{E) and D:l{E), it is useful to 
write the T-matrix for the a-a scattering in the total angular momentum J 

state as 

' sin5^r\E)cos5^r'\E) + cos5^r\E)sin5f'\E) . 

(3.21) 

The parameter "r" in the Coulomb functions Fj[rk) and Gj{rk) is also 
expected to be energy dependent. The form of its energy dependence is not 
easy to determine from known general properties of the system. So we shall 
take the following trial form 

r(E) = roil-zE^), (3.22) 

with "z" a parameter to be determined. It leads to a good fit to the experi- 
mental data. 
J^ = 0"*" states: 

The J^ = 0"^ states are shown in Fig. 3. The ground state of ^Be* is in 
the J^ = 0~^ series and is rather narrow. The next two states in this series 
are broad and have resonant energies at ii^ = 20.2MeV and E = 27A9MeV 
respectively. Since these states are well separated from the ground state, it 
is justifiable to include only the ground state in the R-matrix at low energies 
(E^™ < 12 MeV). 

The best fit (least square fitting) to the experimental a-a phase shifts 
[ ^ in the J^ = 0'^ state is achieved by using the set of parameters shown 
in Table ^. The result of the fitting procedure is shown in Fig. 6. The value 
of the J^ = O"*" phase shift at ii^ = OAMeV given in Ref. [p^ is interpreted 



as 180° ± 0.5 instead of 0° ± 0.5 in order to take into account the fact that 
the O"*" resonance is at about E = 0.2MeV, which is smaller than OAMeV. 



22 

J^ = 2+ states: 

The J^ = 2'^ states are our prime interest. Fig. 3 shows a series of 
experimental low lying states with J^ = 2+. In the energy region considered 
here, we shall assume that only the lowest three states are important, namely 
the states with energy E = 3.04 MeV, 16.63 MeV and 16.92 MeV. Higher 
states are assumed to have negligible influence for ~ 12MeV. The state at 
E = 3.04 MeV is an isosinglet state. The nearly degenerate doublet states 
at an energy near 16 MeV are mixtures of isosinglet and isotriplet. The 
latter component does not contribute to the a-a scattering to lowest order 
in the fine structure constant a = e^ /he. The mixing angle between the 
isosinglet and the isotriplet is nearly 45'' PO, pil]. This mixing angle will 



be determined by fitting the a particle spectra in the /5-decay processes. 
They are not relevant to the phase shift analysis^. In addition, the widths 
of these two states have very small influence on the a-a phase shift in the 
energy region considered in Fig. 7. 



The best least square fit to the experimental a-a phase shifts ||2^ in the 
J^ = 2+ state is achieved by using the set of parameters shown in Table ^ in 
which the widths of the two 16 MeV doublet states are not determined from 
the a-a scattering phase shift but rather from the a-spectra in the /3-decay 
processes. The result of the fitting is shown in Fig. 7. 

It can be noted that the value of r obtained in the J^ = 2+ channel 
is smaller than two times the a particle size. This feature should not be 
regarded as a shortcoming of the procedure. As is pointed out in section 3. a, 
the value of r is not necessarily restricted to be outside of the region where 
the a-a overlap is unimportant. In our approach, r represents the location 
of a zero of the radial wave function before the decay potential AV is turned 
on. From this point of view, there are many equivalent r's. 
JP = 4+ state: 

The J^ = 4"*" series of states are also shown in Fig. 3. The lowest state 



■^■^ This conclusion follows because the a-a phase shift expressed by the R-matrix con- 
tains no information about the resonance other than its width and the location. Without 
further assumptions, the isospin mixing of the resonant states is not easily determined by 
the phase shift analysis in the present approach. 



23 



is at 11.4 MeV. Higher states in the J^ = 4"'' series will be neglected in the 
following. 

The best least square fit to the experimental |^ a-a phase shifts for 
J^ = 4"'' is achieved by using the set of parameters shown in Table |^. The 
result of the fitting procedure is shown in Fig. 8. 

c) The a particle spectra in the /3— decay processes of 
the A=8 system 



The unpolarized differential /5-decay rate follows from Eqs. |2.33| and |2.34| , 
namely 

dW 2 



dEr Svr 



2 



Gj,cos'ec\gArM^JM^Er\A2{Er)\'D{Er), (3.23) 



Eo-Er 



D{Er) = / deeeeJtl - m; 



(Eo -Er- eef F{Z, ee)RT {Er, ee), (3.24) 

FiZ,e.) = J^^, (3.25) 

V - T^P=, (3.26) 

with Z = A and A2{Er) the energy dependent matrix element obtained from 
the shell model amplitude Eq. p.35| using the R-matrix formalism (see the 
following) . 

From Eq. |3.1| , it follows that the matrix element between the initial ^Li 
or '^B state and the final resonant state ^Be* that eventually decays to two a 
particles can be related to a set of shell model matrix elements between the 
same initial state and a discrete set of shell model ^Be* states through 

{Jf\\Oj\\J,) = e-^-/(^^V , Ti^{Er){Jf,n\\Oj\\J,) ^^^^^^^ 

^[Er - En' - Di\Er)f + fj {Erf 

where (J/ ||0j|| J,) denotes any reduced matrix elements, and the index "n" 
enumerates different shell model states with the same parity and spin (total 
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angular inomentuin) J/. The common phase factor for each of the terms in 
the summation is written as a multiphcative factor of the whole expression 
on the r.h.s. of Eq. p.27| with the vertex function T^ (E) real. This ex- 
pression is, in principle, exact provided a complete set of shell model states 
is used and the energy dependent vertex function r;((i?), level shift D:^{E) 
and half width li{E) are calculated from Eqs. p.6| , p.7| and p.8| . Due to 
the lack of knowledge of the decay potential AV^ and the complexity of the 
problem, a microscopic calculation of V'^{E), D^{E) and li{E) is beyond 
the scope of this work. Instead, we shall use the same parameterization of 
the R-matrix as in the a-a scattering case. Consistent with the treatment 
of the a-a scattering phase shift, only the 2"*" shell model states at E=3.55 
MeV, E=16.19 MeV and E=13.37 MeV are used to evaluate the shell model 
transition amplitude. The energies of these states are taken from the fit of 
the a-a phase shifts. (This is reasonable, since it is expected that there are 
level shifts for the shell model states after switching on the decay potential 
A\^). The shell model state lying at E=17.45 MeV is not included in the sum 
of Eq. |3.27 . Two additional parameters in the /9-decay processes that do not 



appear in the description of a-a scattering processes are needed. The first 
parameter is the isospin mixing angle between the isospin doublet members 
at £" ~ 16 MeV. We shall define the mixing angle as 

1^ = 16.63) = cos0|T = l) + sm0|T = O), (3.28) 

|E = 16.92) = -sm0|T = l) + cos0|r = O), (3.29) 

where | T = 0) and | T = 1) denote the isosinglet and isotriplet shell model 
states respectively. Since only the T=0 components within each state of the 
doublet can decay into two a particles, the contribution of the T=l part of 
each state is not included in the R-matrix. In addition, there are relative 
signs between 1^=2'' (3.04), 1^=2'' (16.63) and 1-^=2^(16.92) that cannot be 
fixed from the sum rule relation Eqs. ^^ and |3.5| . They have to be determined 
in the fitting. The second parameter is related to the inclusion of the potential 
scattering in the R-matrix formalism for the /3-decay processes consistent 
with the similar development for the a-a scattering. This parameter is 
necessary, since in the fitting of the a-a scattering phase shift, it is natural to 
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separate the potential scattering between the two a particles at long distances 
and the confining strong and Coulomb resonant interaction at short distances 
(characterized by the parameter r). Similar to Eq. |3.21| we shall write the 



T-matrix in the presence of the potential scattering between the a particles 

as 



^{pot)(£;^ 



(pot) { rp\rr { rp\ I „„-^x(P°*)/ 



T{E) = e''f ^''>[cos6)r'{E)To{E) + stn6T'{E)To{E)\. (3.30) 

Tq{E) has the form of Eq. |3.1| but with the full vertex function r:^{E), the 
level shift function D:^{E), and the half width function I^iE) replaced by the 
ones obtained from the phase shift fitting procedure outlined in the previous 
section. Tq{E) is given by 

UE) ^ Ci:g|| ^:g:g|g (";eH^I^.>. (3.31) 

with constant C a parameter to be determined. The magnitude of C is a 
rough measure for the direct weak transition from the initial ^Li or ^B to the 
two a scattering states. 

The two new parameters C and 0, along with the widths of the isodoublet 
states are adjusted so that a fit to the a particle spectra for the (3^ / (3~ decay 
in A = 8 system is achieved. The resulting values for C and are 

(/> = 51.8°, (3.32) 

C = 5.7 X IQ-l (3.33) 

The widths of the doublet state are 

w;o(16.63) = 0.206MeV, (3.34) 

m;o(16.92) = 0.303MeV. (3.35) 

In addition, relative to 1-^=2'" (3.04), the signs of r-^=2+ (16.63) and r-^=2'" (16.92) 
are negative. Before ending this subsection, let's compare the isospin mixing 



angle obtained here with the one obtained in Ref. [31|, namely = 50.3°. 



Thus, the value for the isospin mixing angle obtained here agrees well with 



that of Ref. pl| 
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The results of the theoretical calculation, with the nucleon weak form 
factors given in Table 7, are compared to the experimental [^ data in Figs. 



9 and 10. The dashed curve represents a calculation without the final state 
potential scattering, namely 5^^°^'' = and the solid curve is the full calcu- 
lation. Figs. 11 and 12 are linear scale curves corresponding to Figs. 9 and 
10. 

In obtaining the excellent fit to the a particle spectra shown in Figs. 9 
and 10, an upward energy shift of 0.17 MeV has been made to the theoretical 
expression Eq. p.23| . It is as if only the lowest resonant energy J^ = 2"*" state 
were shifted upward by roughly 0.17 MeV. This shift spoils the agreement 
between theory and the experiment both in the shape of the a spectra and 
the phase shift. Furthermore it is much too large to be attributable to an 



experimental energy calibration ^3[. We do not understand its origin, only 
its necessity. In Warburton's fit to the a-a scattering phase shift and the a 
particle spectra of the /3-decay processes of the A=8 system, a shift in the 
resonant energy of order 0.07 MeV was also required 0]. Henceforth we 
shall use the R-matrix development without this shift. 

It may be worth mentioning that the Barker and Warburton R-matrix 
treatment, which is different from the one used here, obtains an equally good 
(if not better) fit to the a-spectra EM without the need of low lying "in- 



truder" states if the value of Re is chosen to be around 4.5 fm. Their R-matrix 
formalism does not, however, allow a straightforward and consistent exten- 
sion to the forbidden transition matrix elements. On the other hand, in the 
R-matrix developed here, the allowed and forbidden transition amplitudes 
are treated on an equal footing, so that the study of PTRIV phenomenae, 
which are intrinsically due to forbidden processes, can be studied with the R- 
matrix determined in the previous sections. This kind of "predictive power" 
is the main reason that we adopted the present approach. 

4 PTRIV due to the final state interactions 
between two final a particles 
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a) Energy dependence of the PTRIV response func- 
tions 



Eqs. |2.48H2.51| show that the dependence of i?f ^-^^^ and R^'^^^^ on the en- 



ergy of the charged lepton energy eg is quadratic. However, their dependence 
on the two a particle relative energy is not trivial. First there are strong E^. 
dependence of ?76(2+, 0+), 7/9(2+, 0+), 7/6(2+, 4+) and 7/9(2+, 4+) . Figs. 13-16 
give the energy dependence of these quantities. Functions 775(2+, 0+) and 
7/9(2+, 0+) are dominated by the J^ = 0+ resonance located at ii^ = 0; they 
have a zero near E' ~ 3 MeV. 7^6(2+, 4+) and 7/9(2+, 4+) are dominated by the 
J^ = 4+ resonance located at i? ~ 11 MeV; they do not possess any zero in 
the energy range 2-5 MeV. Also, they depend on factors sin{62 — 60) and 
sin{62 — 64), which originate from the phase factors in the wave functions 
(see Eqs. ^]6| and p.3(]| ). There is an energy near E ^ 3 MeV where the 



J^ = 0+ phase shift, 6q, is equal to the the J^ = 2+ state phase shift as can 
be seen by a comparison of Fig. 6 and Fig. 7. Therefore, there is a zero in 
sin{62 — Sq) . However, sm((52 — 54), does not have a zero in this energy range. 
In a test of genuine time reversal invariance, which is associated with the 
time reversal invariance of the underlying Lagrangian, experiments should 
be designed to select events in which the PTRIV contributions are small. In 
Figs. 17 and 18 the value of Rf^^^^ and i?f-^^^^ are plotted against Er and 
eg for the /3+ decay of A=8 system. It can be seen that both of R^^^^^ and 
jlPTRiv pa^gg through zero for several combinations of Ej, and eg , along a line 
located at Er ^ 3 MeV. R^^^^^ is small at low e^ energies (eg) and increases 
with ee; it changes more quickly with E^. R2^^^^ changes less quickly at 
eg ~ 5.5 MeV than at other values of eg and possesses zeros somewhere near 
Er ^ 3 MeV and eg < 8 MeV. Figs. 19 and 20 are 3-dimensional plots of 
the values of R^^^^^ and i?^-^^^^ against Er and eg for the /3~ decay in the 
A=8 system. The general features of R^^^^^ are the same as those for the 
/?+ decay, i?^-^-^^^ is somewhat different from the corresponding /3+ decay 
due to the v — a interference effects. There still are zeros around Er ^ 3 MeV 
at different ee, but R^^^^^ turns negative at high values of Er- It changes 
slower with Er when eg ~ 5.5 MeV compared with its value at other values 
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of ee in the energy region (of Er) considered. 

b) PTRIV in case the charged lepton energy is not 
measured 

In case the charged lepton energy is not measured, an average over its energy 
distribution has to be made. From Eq. |3.24| , it is natural to define the 
averaged PTRIV response function as 

iV Jnie 

(4.1) 
N = / de,eJel-ml{Eo-Er-e,YFiEr.,e,), (4.2) 



with i = 1,2. 

The Er dependence of the average PTRIV response functions are shown 
in Figs. 21 and 22. 

5 e^— nucleus Coulomb interaction and related 
PTRIV 

The treatment of the Coulomb effects in the /3-decay processes has a long 
history I^S], ^ ^ . However the coordinate space version of the treatments 
are not easily adapted to the formalism developed in this study. We shall 
use a somewhat different approach. 

Although it is not necessary to consider the full effects of the charged 
lepton-nucleus scattering for the purpose of this paper, we nevertheless treat 
them here. The particular aspect of the Coulomb scattering effects of the 
charged lepton studied in detail here is the Coulomb scattering induced 
PTRIV effects in the /3-decay processes of the A=8 system. Similar effects 
were studied by Jackson, Treiman and Wyld for general allowed /3-decay pro- 



cesses 38 . The Coulomb PTRIV effects were later related to e -Coulomb 
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phase shifts |^, similar to the relationships developed in the previous sec- 
tion. When we go beyond the non-relativistic limit, infinite series of partial 
waves are required (see later sections of this chapter) . For our specific prob- 
lem, we shall not start with the spherical basis to calculate the contribution 
of the e^-nucleus Coulomb scattering to the PTRIV, but rather make use of 
the integral form of the Dirac equation for the charged leptons. These sets of 
integral equations are relativistic Lippmann-Schwinger equations, with the 
boundary conditions built-in from the beginning. 

a) Modified weak interaction Hamiltonian due to the 
e^— nucleus Coulomb scattering 

When the Coulomb scattering of the e^ due to the presence of the charged 
nucleus is considered, q can have a distribution of values associated with 
different kg at the point of the e^ emission. If we take into account the 
Coulomb scattering, the weak Hamiltonian can be written as 

Hw = ^/^(ke,k.|j.(q)|0)J'^(-q), (5.1) 

where the helicity indices for the leptons have been suppressed. 

The matrix elements of the leptonic weak current operator in the /3-decay 
processes can be written as 

(k,;k,|j^(q)|0) = Jd'xe''^--{K;K\j^i^)\0) 



^plr\{Kh,{l - ^')v,iK) for p- decay 
MKh,{l - 7^)^Sk(ke) for /3+ decay 



(5.2) 



with Mi/(k) and Vi/(k) massless neutrino and antineutrino spinors respectively. 
Eqs. |5.1| and |5]^ can be compared to earlier work on the Coulomb corrections 
to the allowed /3-decay in the elementary particle approach ||D, |T, ^. The 
Fourier transformation of the e^ wave functions are given by 

^|Je-,l.M) = /rf3a:e-*'^^-^Ve-,k.(x), (5.3) 
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V'e+.k.lk'J = y'ci3a;e*'''^->e+,k.(x). (5.4) 

b) e* wave functions and their relativistic Lippmann— 
Schwinger equations 

The wave functions for the e^ obey Dirac equations for relativistic charged 
fermions in Coulomb fields 

{ho + V + Pme)iJe- = ee^e-, (5.5) 

{ho - V + Pme)i'e+ = ee^e+, (5.6) 

where V = — Za/|x| in coordinate space and 

ho = —id ■ V + /5me, (5.7) 

with a* = 7^7* and /3 = 7*^. ho is the free Dirac particle Hamiltonian. A 



very detailed study of the scattering state solutions to Eq. ^]5| or ^]6| can be 
found in the book by Greiner, Miiller, and Rafelski |^ . 



For our purpose, it is more convenient to use a relativistic Lippmann- 
Schwinger equation for the e^ wave functions, namely 



^^^^^(x) = M(ke)e*'^-^+ /t/V(x| i^Ti^)|x>(ke)e^''-^', 

(5.8) 
V'SkW = ^(ke)e-'^--+/rfV(x| l^_T(f|x>(k,)e-^'^-''', 

' J — 6(, — llo i ^C 

(5.9) 

where eg > 0. The free Dirac Hamiltonian ho and the matrix elements (x | 
. . . I x') are 4x4 matrices in Dirac space. In fact, Eq. |5.9| corresponds to the 
negative energy solution of Eq. |5.8| , but in order to use field theoretical lan- 
guage, ee is restricted to be positive and equations for e~ and e"*" are written 



separately. Similar to the free Dirac spinor [^, the positron wave function 
can be obtained from the electron wave function by a charge conjugation and 
a change of sign of the Coulomb potential, namely 

V'SljK-.] = C7VI?* [v;+], (5.10) 
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where K+ is the Coulomb potential experienced by a positron and [V] de- 
notes the potential used in solving the corresponding wave function. Due to 
the long range nature of the Coulomb potential, the Lippmann-Schwinger 
equation may not be well defined. This difficulty is overcome by introducing a 
screened Coulomb field, which is the case in any realistic experimental situa- 
tion. The "pure" Coulomb case is realized in the limit of an infinite screening 
length. In most of the following formal discussions, we shall assume that the 
above limiting procedure is taken. 

The Dirac spinors M(ke) and t'(ke) represent a free electron and positron 
respectively; they satisfy 



(^-me)M(k) = 0, 



(5.11) 
(5.12) 



T'„_ /„ , is the T-matrix of the e^ in the Coulomb field of the nucleus under 



" e^/e 



.(±) 



(±) 



study. T _ and TV satisfy equations 



rp{±) 



T(±) 



v + v 



I 



Ce - /lo ± ie 



i-v) + {-v)- 



1 



t{±) 



(5.13) 
(5.14) 



-Ce ~ ho±ie 

In momentum space, the Lippmann-Schwinger equations for an electron and 
a positron in a Coulomb potential become 

1 



C(k'e; 



^SUk: 



{2nf6{K - K 



{2nf5iK - K. 



a ■ kg — PrUe ± ie 



1 



— Cg + a -k' — Prrie ± ie 



Ti"^(k;,k, 



Te+ (^^6 5 kg 



Define the pair of operators 



<±) 



Ol?(k'e,kg 



ee — a ■ kg — j3me ± ie 
1 



TgV"^(k'e,ke), 






M(kg), 
(5.15) 

(5.16) 



(5.17) 
(5.18) 
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The conjugate operators are 

o«(k.,ki) ^ ri?(k.,k;)— -^;^L__, (5.19, 

oS'iKX) = ^SHkk)^^^^^^^^^^—:^. (5.20) 

These operators will be used when the matrix elements of the leptonic weak 
current operators are studied. 



c) The leptonic tensor 



Using operators given by Eqs. |5.17 - ^^ the leptonic tensor corresponding 



to the trace of the bilinear products of the leptonic weak current operators 
becomes 

Tr/(q)/^(q) = (27r)65(q + ke + k,)5(q' + ke + k,)L'^'^ 

+ (27r)35(q + ke + K)A^'{c() + (27r)35(q' + k^ + k,)5'^'^(q) 
+C^'^(q,q'), (5.21) 



where L^"" is given by Eq. |J. A^"' , B^" and C^"' are 
For the P~ decay 

A'^'^iq') = 2TrOe-(-k,-q',ke)(|^e + me)7^|£,7^(l-7'), (5.22) 

B'^'^iq) = 2Tr{l(:, + m,)0,-{K,-K-ci)nt,Y{l-l'), (5.23) 

C^'^(q,q') = 2TrO,-(-k,-q,ke)(|^e + me)Oe-(ke,-k,-q)7%7''(l-7^). 

(5.24) 

For the j3~^ decay 

A'^^(q') = 2Tr(|£e-m,)a+(ke,-k,-q')7%7'^(l-7'), (5.25) 

B^'^iq) = 2TrO,+ {-K-ci,K){lte-m,)Yhr{^-7'), (5.26) 

C^'^(q,q') = 2TrOe+(-k.-q,ke)(|^e-me)Oe+(ke,-k,-q')7%7^(l-7'). 

(5.27) 



33 



d) On— Shell Coulomb T— matrix to all orders in the 
fine structure constant a 

Rigorous expressions for A^^ ^. . . ,C^^ involving the T-matrices (on-shell and 
off-shell) of the charged lepton in the hadronic Coulomb field are expected 
to be very complicated. Due to the smallness of the fine structure con- 
stant a = 1/137, it is reasonable to adopt an approximation in which only 
terms to the first order in Za are considered. This however can not be done 
straightforwardly by truncating the Lippmann-Schwinger equation, since the 
truncation leads to infrared divergences for many of the PTRIV observables. 
Even when a screened Coulomb potential is used, the value of the T-matrix 
in the forward direction cannot be reliably obtained using the lowest order 
one, which is real, when the screening length is much larger than the hadronic 
scale. In fact, an all order (in a) calculation yields a rapidly varying phase for 
the T-matrix in the forward direction that cancels the infrared divergences 
in the physical observables. This phase can be calculated analytically in the 
non-relativistic limit ||2^ (i.e., when the momentum of the charged lepton is 
much less than its mass). Away from the non-relativistic limit, the Coulomb 
T-matrix does not have a closed form. It is however possible to express it 
in terms of the Coulomb phase shifts obtained in Refs. ^, ^ and more 
recently in e.g. Ref. [^. The first order approximation in Za is made only 



at the final step of the calculation. For the special case of the /3-decay pro- 
cesses of the A=8 system, the electron (positron) energy is much large than 
its rest mass. This allows us to use the T-matrix in its extreme relativistic 
limit, where helicity is conserved. In this limit, it is easy to show that the 
Dirac structure of the T-matrix is very simple. It can be written as 

t = t2®l, (5.28) 

with "1" an unit 4x4 matrix. The matrix elements of t between states with 



(5.29) 

(5.30) 
(5.31) 



definite helicities have the form 




t++(k',k) 


o 

= cos-t2{e,e, 


r-(k',k) 


= t++(k',k). 


t+-(k',k) 


~ r+(k',k) = 
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where "+" represents helicity "up" and "-" represents helicity "down" and 
{9,(f)) are the polar angles between k and k'. Quantity t2{e,9,(j)), which is 
independent of 0, is related to the Coulomb phase shifts through 



F Ik 



1 oo 

(5.32) 



with 5±^ Coulomb phase shifts of the doublet states with total angular mo- 
mentum j = |k| — I and with opposite parities, -P|(x) = dPi/dx and Pl{x) is 
the Legendre polynomial. Due to the helicity conservation of the T-matrix, 
the phase shifts of states with the same total angular momentum and with 
opposite parities are the same. Using the fact that 6^ = 5-ki Eq- |5.32| can 
be written as 



1 oo 

t2(ee, 9) = -^^ Y. [e^^'^smd^ - e^''^~'sin5^.i)\ P'^_^{cos9). 

(5.33) 



f Ik 



PTRIV correlation observables depend on the moments of the real part of 
the on shell Coulomb T-matrix if a first order approximation in Za is taken 
in the final results. These moments are defined as 



Ik P /•! 



Zaep 



I dxx''Re[t2{ee,x)], (5.34) 



with X = cos9. Taking a first order approximation in Za for the Coulomb 
phase shifts 6^, an can be expressed as 



n „i 



^(l + (_l)-+-)A/,-n^ / dxx-~'P^^,{x)Af,, 

K=l K=l"'-1 

(5.35) 



A/.^i = 


1 1 r 
2k{k-1) ' 2(fi:-l)! ^^'^ 


-in^- 


-2)- 


- k{k - 


- 1) + 2e(2 - K)\ , 
(5.36) 


A/i = 


1 








(5.37) 
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where e = 0.57721566 ... is the Euler-Mascheroni constant and the Coulomb 
phase shift obtained in Ref. 



has been used. The first few useful moments 



are 



ao 
ai 

as 



.(+) 



.(-) 



with 



c(+) + 1, 



2 E ^/- 

K=2,4,... 

oo 

c- = 2 E AU 

K=l,3,... 



.(+) 



.(-) 



(5.38) 
(5.39) 
(5.40) 
(5.41) 

(5.42) 
(5.43) 



The numerical value of c'-^-* and c^~^ are given by c'-^-* = —9.942 x 10~^ and 
c^^^ = 9.942 X 10^^. It can be easily seen that all these moments diverge 
if a first order (in Za) approximation for the Coulomb T-matrix of the 
e^ is usedp^ due to the singularity of the first order T at x = 1. The 
results of the Coulomb T-matrix for the e^ to all order in Za suppresses 
the divergence at x = 1 and results in finite moments a„. In addition, it 
can be shown that, except for some special combination of the moments, a„ 
contains contributions from states with all possible total angular momenta 
when the mass of the e^ is not much larger than its 3-momentum (i.e. away 
from the non-relativistic limit). However, in the extreme relativistic limit, 
the average effects of these states, represented by c^~^^ for even n and c'-"^ for 
odd n, are small and the sum of c'-"*"-* and c^~^ seems to be vanishingly small 
in a numerical evaluation of the corresponding series. 



^■""^Some combinations of the moments remains finite even in the first order approximation 
for the Coulomb T-matrix. 
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e) Dispersive part of the Coulomb corrected leptonic 
weak current 



Only the dispersive parts of Eqs. |5.22| - |5.27| will be studied in this paper, since 
they alone contribute to the PTRIV signals. In addition, to the first order 
in Za, the dispersive parts of Eqs. p.22| - |5.27| originate from the imaginary 
part of the free Dirac particle propagators 



1 



P 



' ho ±ie 
1 



ho 



P 



-Ce -ho±it -Ce 



h 



=Fi7r(5(ee-/io), 
=F Z7r(5(-ee - ho 



(5.44) 
(5.45) 



where "P" denotes principal value. 
For /?~ decay 



D±\disp{.^.,^) 



and 

for (5^ decay 

Df|,,,,(q,q') 



vr 



|ke|/?e[tr(ke, KWrihrK T Krhh'hYii - 1 



(5.46) 



(5.47) 



where D'^ is related to A^'^ and B'^'^ given in Eqs. p. 22 and 5.23 or Eqs. 
5.25| and ^.26| through 



^f(q,q') = j^Jdhkl[A^''{q)±B^''ic^)]6ie,-e,),, (5.48) 



with q = -kg - k^, q' = -k"e - K, K = \k'J = |k"e| and e^ = ^k^ + m^. 

The symmetric (with respect to A'^'^ and B'^'^) and dispersive parts of 
the final e^-nucleus Coulomb scattering contribution to the leptonic tensor 
given by -D+'^(q, q') has opposite time reversal transformation properties to 



37 



the non-interacting part of L^'^ given in Table 4. Therefore the T-odd kine- 
matic functions of rank 2 are given in Table ^, which can be compared to 
the ones given in Table |I]. The antisymmetric (with respect to A^'^ and B^^") 
and dispersive parts of the final e^-nucleus Coulomb scattering contribu- 
tion to the leptonic tensor given by D^'^(q, q') has the same time reversal 
transformation properties as the non-interacting part of L^'^. Both of these 
kinematic functions contain an extra variable kg which has to be integrated 
out. 



f) PTRIV correlation observables 

The dispersive contribution of the e^-nucleus Coulomb scattering amplitude 
to the T-odd differential /3-decay rate dW of the A=8 system can therefore 
be written as 

,1 



dW = 2n5{E,-Ef- 

d^q d^q' 
(27r)3 (27r)3 



Cp - e,. 



-Gicos'e, 



Tr 



jMt(_q)jv(_q')lTr[jt(q)j,(q' 



disp 



(5.49) 



By using the hadronic trace formula [0, and after some algebra and inte- 
gration over dQ]^/, it is found that in addition to Eqs. p.46| - |2.5l| , there are 
other terms that contribute to dW^^^^^ . Both D':^ and Z^^'^ contributes to 
the PTRIV correlation observables. The PTRIV observable in the /3-decay 
processes of the A=8 system can be written as 

dw^TRiv ^ dWoT2{k): 

(rPTRIV I nPTRIV\ 



+ RZr)^eky<k + {R 



PTRIV , tdPTRIV\ C t ^C 



em, 2 J 



(5.50) 



where the tensorial contraction between the second rank Cartesian tensor T2 
and two vectors by symbol ":" is defined by Eqs. |2.36| and p.37| . It is found 
that to order k^ in the hadronic response function are 



R 



PTRIV 
em.l 



TZa 



^A ±A + /„ + /. + i 




/iw + 
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W^'~ 



ee(A - 6e 
M2 



-M X 



1 /Is, V2 1/10 2 1 

-30 VY^^^^ ^ W^^ ^ 60 VY^^ ^ 1^7^^^° 



(5.5i: 



R. 



PTRIV 
em, 2 



TZa 



ee(A-ee)^ / 12 fs" ^2 1 /lO 4 1 



..^^^4;;^.3-^^..)x 




72 1/10 2 1 

30 V 1^'"^' ^ ^"^^ ^ 60 VT''' ^ 1577'''° 



where the reduced response functions rji are given in section 2, and 



(5.52) 



/?2 = 3a3 + 3a2 — ai — oq 
/^a = 3a2 + 2ai — oq = 4 — 2e 
/?4 = 2ai + 2ao = 1 - 2e. 



14 



(5.53) 

(5.54) 

(5.55) 
(5.56) 



In order to test whether or not the underlying Lagrangian of the system is 
invariant under time reversal, one should concentrate on events in which the 
PTRIV contribution is small, so that the genuine TRIV can be extracted. 
Figs. 23 and 25 are 3-dimensional plots of i?^f^^ and Rem2^^ fo^ the 
/S"*" decay process of the A=8 system. R^^^^^ for the /J"*" decay process, 
which contains the first forbidden contributions (namely contributions of 
order k/M) obtained by Holstein [Q, is of order 1 — 6 x 10^^. Rem 2^^ 
for the Z?"*" decay process, which contains only second forbidden retardation 
terms (namely, terms of order kR with R some typical radius of a nucleus) , 
is of order or less than 5 x 10~^. Figs. 24 and 26 are 3-dimensional plots of 
RZ^^^ and RZ^^^ for the /?+ decay process of the A=8 system. R^l^^^ for 
the f3~ decay process is smaller than the corresponding one in the /3+ decay 
process due to the destructive interference of fc and /m. It is of order or 
less than 10~^. -R^f^^ for the /3~ decay process depends only on the second 
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forbidden retardation terms. It is of order or less than 2 x 10^'^. From these 
graphs, it can be seen that the observables proportional to kj^ke x k^ in both 
f3~^ and f3~ decay processes contain smaller leptonic Coulomb interaction 
induced PTRIV than other observables; thus it is advantageous to use these 
observables to study the question of genuine time reversal violation. The Er 
dependence of R^^^f^^ (i=l,2) is determined by the energy dependence of r^g, 
r]s,- ■ ■ ,f]io- As we can see, the PTRIV signal due to the e^-nucleus Coulomb 
interaction reaches its maximum in the energy region Er ^ 3 MeV, which 
is in contrast to the hadronic PTRIV. Therefore, although the e^-nucleus 
Coulomb interaction is suppressed by a factor of Za, the magnitude of the 
PTRIV effects resulting from the hadronic strong interaction are smaller than 
those due to the e^-nucleus Coulomb interaction in the energy region Er = 3 
MeV. Of course, the hadronic strong (and Coulomb) PTRIV can be bigger if 
one goes away from the the Er = 3 MeV region, in which different hadronic 
PTRIV amplitudes cancel each other. 

6 Conclusions and Acknowledgments 

In this paper, we have examined a time reversal test in the /9-decays of the 
mass 8 nuclei. An R-matrix treatment was used, with parameters related to 
the a-a scattering phase shifts and allowed /9-decay rates. The e^-nucleus 
Coulomb final state interactions are included. False time-reversal violation 
(PTRIV) signals could arise due to these e^-nucleus Coulomb as well as 
the strong interaction effects in the A=8 system. The strong interaction 
induced PTRIV is much smaller than the e^-nucleus Coulomb interaction 
ones because the former is due to several second forbidden contributing terms 
that nearly cancel each other in the Er ~ 3-5 MeV region. In the energy 
region close to the lowest 2+ state around 3 MeV these contributions can 
be further minimized because they have zeros. The e^-nucleus final state 
interaction effects are < 7 x 10~^. These effects have not been seen so far. 

The authors are particularly grateful to Dr. Ludwig De Braeckeleer for 
many useful discussions of the experiment and associated theory. We also 
thank Eric Adelberger for helpful comments. 
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Appendix A: Hadronic Response Function in 
Terms of the Reduced Matrix Elements of Mul- 
tipole Operators 



Following Ref. [18 
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(A.l) 
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Jf 


J'f 
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Ji 
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>F„,{JJ';apT), 


(A.3) 


J 


J' 


o 







where x = \/2x + 1, in order to simplify the notation, a sum over indices not 
on the l.h.s. of the above equations is implied. It is found that 



F„(JJ']a) 






(A.4) 
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F±(J/;ap) = fsviJJ';ap)±fUJJ';(^p), (A.5) 

1 /' J' 

r p a 






and 






J' I p 
J 



(A.6) 

(A.7) 
a/C}, (A.8) 



and 



«j-i 



^j+i 



-^y^J^ (vOn^j + VJCj) , (A.9) 

-lYMj, (A.IO) 

-l)'^j^ (v^^J - ^^JTTCj) . (A.ll) 



The reduced matrix elements Cj, Cj, Sj and M.j oi the hadronic charged 
weak currents can be decomposed into 



Cj 


= Cj + C]={Jf \Cj\ J,), 


(A.12) 


Cj 


= Lj + L''j={Jf\\Lj J,), 


(A.13) 


Sj 


= Ej + E] = {Jf\\ff\\J>), 


(A.14) 


Mj 


= Mj+M'j = {jf fr'\\J^), 


(A.15) 



where superscript "5" indicates that the corresponding matrix element is 
originated from the axial vector current and the rest of the matrix elements 
are originated from the vector current. 
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Appendix B: Power Expansion of the Hadronic 
Response Functions for A = 8 System in Terms 
of Momentum Transfer 

The relation between Ri . . . R12 and Rgs ■ ■ ■ Rw can be found. For T-even 
observables, only the final 2"*" state contributions are considered. While for 
T-odd observables, we consider contributions from the final 2+ state, and 
from the interference between the final 2+ state with the final O"*" and 4"*" 
states. In the following discussion, we shall introduce set of response func- 
tions Rss, ■ ■ ■ , Rvv that are related to Rgs, ■ ■ ■ , Rw through 

p 

^* = I — \2\A 12 ' {i = ss,...,vv) (B.l) 

where A2 is the Gamow-Teller matrix element defined in Table ^ With 
these set of response functions, the Cartesian ones can be found. The mass 
of the charged leptons rUe can be neglected; we do so in the following. Those 
correspond to cr = and T-even observables are 

4°^ = 5RJ0) + 5— i?i+)(Ql) - ^VSR^JOOO) ± 5^2 ^' ~ ^^ RwiOll) 

+ fyf^i?..(022), (B.2) 

i?f = 5Rss{0) + 5-i?i+)(01) + ^i?„„(000) T 5V2^^^^Rwi0n) 

+10^^^^^i?..(022), (B.3) 

Rf^ = 5^-'-^RU022). (B.4) 

Those correspond to cr = 2, T-even, which contain only the final 2+ state 
contributions, are 

i?P^ = -^4^..(2) + ^^i?W(21)T5./^^i?L)(22)±15j^i?„.(211) 
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+ ?\/44^-(220) + ^v^^i?..(222), 
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-%'-^Rss{2) - 25 J;^^i?..(220) - ^v/30^i?„„(222), 
5 



7' ""/t2 
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Tl5W-^i?..(211) + A^W^^i?,,(220) + -^v^^i?„„(222), (B.7) 



(B.8) 
(B.9) 



/30^i?,,(222). 

Those correspond to o" = 2, T-odd, which contain only the final 2"*" state 
contributions are 
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15 e 
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Those correspond to a = 2, T-odd, which originate from the interference 
contributions between final 2+ state and final O"*" state, are 



(2) 

Rg = sin{62 - So) 



2e.{ee-e.) 



K^ 



Re {Rit\22))^^ ± 15^^i?e (/2(;)(21) 



20 



T 



lO^'A^^^Re (i?,„(221))^^ - W\^-Re (i?.„(212) 
cos{62 - So) 



20 



^ /2 M^„„ (^,,,,22))^^ ± 15 /2!,,,,„ (;,,;,,2i))^^ 



TlOv^ 



eeUe - ^u 



K^ 



Im (r,,{221))^^ - 10i/|-/m (i?.„(212) 



20 



5 K 



20^ 
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ez/(ee - Cj, 



K^ 
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'''t3 
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Those correspond to o" = 2, T-odd, which originate from the interference 
contributions between final 2+ state and final 4+ state, are 



R 



(2) 
11 



sin{62 — 64 



30i 



3eJep, - e 



eV"=e "=i// 



7 ^2 



5 e^ 



i?e (i?W(22))^^ ± 18Wpi?e (i?M(21) 



24 



T30i/^ '^^'% '"^ Re (R..{221))^^ - gJ^'-Rc (i?..(212) 
V7 K^ ^ ^24 V7/t^ 



24 
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The subindices "20" and "24" in Eqs. [RT3| - pl6| denote they are the in- 
terference terms between the 2"'" state and 0"*" or 4"^ respectively. Response 
function R^s, ■ ■ ■ , Rm are related to the reduced matrix elements of the mul- 
tipole operators of the hadronic charged weak currents; they are given in 
Appendix A. Using the definition for the static multipoles given in Tables ^ 
and 1^, the response functions are ready to be expressed in terms of /i, /m, 
It, fc and r/fc (fc = 1, . . . , 10). (See Eqs. 2.41-2.45.) The resulting T-even 
response functions have the following form 



m 



(0) 



p(0) 



^4eep / , ^5 /2A V2A 

2el ^ I, [2 2^2 

4e, 



2V2 A 



^8 



01/ /2^ 2v/2A 



(B.16) 
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V7- 



2v^ A 
Is"!?'"" SVSSM''" ' lOSv^M' 



r% 



2A 

e' o /^2 2v^ 



4 / 2 A 76 A 

:V9 + ——^77^0 



76 



^8 



3V 35 



m 



105^7 



^10 



:b.2i) 



^ / /5 . . ^ f3"A v^A 



2 A 



4 A 



35m''^ + V7m''^° 



ee / . /5 A V2A /2 A 



d(2) 
-K4 



48 A 
'3577M 



\ (A-ee)2 / 1 /7 ^ /2 2 



i? 



(2) 



M^"^" (^W T^^^^^ + T^v y^^ " ^177^^° 

ee /2A 2^2 A 12 A \ 



;B.22) 
;B.23) 



50 



K 



(2) 
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-^8 + 
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Before ending this appendix, some more details are worth mentioning. The 
general form of response functions Eqs. |B.2HB.10| depend on kg ■ k,^ though 
their dependence on k. When the power expansion in terms of k to order 
0{k^) is performed, these dependence are extracted to make the response 
function defined in Eqs. p.l7| - |B!25| independent of kg ■ kj,. Therefore R\ 
(i=l,2,3) and R^^^ (i=l,. . . ,6) defined in Eqs. §]2]^B]10| are slightly different 
from those defined in Eqs. [B.17| - py25| , which are independent of kg ■ kj,. 
However, the differences are small and only depend on 777. 
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Appendix C: Non— relativistic Reduction of the 
One Nucleon Matrix Elements of the Hadronic 
Weak Current Operators 



p(±)(x) 



pW^(x) 



J(±)(. 



F,^^rW«<5(x-r,) 



j=i 



T'^T.r^^H^M^)■V6{^-r^^ 
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+ ||E^^^^(OVxa(zMx-r,), (C.3) 

^^^Hx) = (^A±^^T)ErW(0^(0^(x-r.), (C.4) 

where A = Ei — Ef is the difference of the initial and the final energies of 
the hadronic system, G\f = Fj^ + 1 = 4.71, terms that depend on gp are not 
included since their contribution to the differential decay rate is extremely 
small when the charged lepton polarizations are not detected 



Appendix D: R— matrix theory for the (3-deca.y 
processes of the A = 8 system 

We only present an outline of the theory here; for more details we refer the 
reader to Refs. [19,26]. 

To the first order in the weak interaction Hamiltonian H^ and all orders 
in the strong and electromagnetic interaction Hamiltonian Hg and -ffem? the 
S'-matrix for the ^B or ^Li -^ 2a transition can be written as 



oo 



Sfi = -i dt{2a;eiy\U{oo,t)H^U{t,-oo)\(f)i), (D.l) 

J —oo 

where f/(ti,t2) is the full propagator of the system under strong and elec- 
tromagnetic interactions and | (pi) is either ^B or ^Li. With a proper phase 
convention for | (pi), which is regarded as an eigenstate oi H = T + Hs + Hem 
(T is the kinetic energy operator), 

/oo 
dte-'^^\2a; eu \U{oo, t)Hj (pi). (D.2) 

-oo 

A complete set of states can be inserted between the propagator and the 
weak interaction Hamiltonian H^j, 

/oo 
rfte-*^»*5](2a;ez/|f/(oo,t)|/)(/|i7^|0i), (D.3) 
-oo J 

which is illustrated in Fig. 1. Since the matrix elements of the weak inter- 
action Hamiltonian Hyj are large only when the hadronic part of the state 
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I /) is comparable in spatial extent to that of | 0j) p^ the summation in Eq. 
p.tj| can be saturated to a satisfactory precision by a subset of states with 
comparable spatial extensions to | 0,) . These sets of states can be generated 
by choosing eigenstates of a Hamiltonian that confine them in a region with 
a size that is comparable to | 0j) . It is natural to choose the direct product of 
shell model states | n) and the leptonic states | ev) . The shell model Hamilto- 
nian is so constructed that the energies of its eigenstates are approximately 
located at the resonant peaksp^ of the ^Be* system. At low energies, only a 
very few such shell model states are important, since the higher excitations 
of the shell model Hamiltonian have energies well separated from the ground 
state. 

The strong and electromagnetic interaction Hamiltonian can be written 
as 

Hs + H,m = Hshell + ^V, (D.4) 

where Hgheii is the shell model Hamiltonian and AV is the residual interaction 
that is responsible for the decay of the shell model ^Be* states. The decay 
of the localized state | n) , of interest in this paper, to an asymptotic two a 
state can be parameterized in terms of the i?-matrix ||2^ as 

R2a,n{t) = {2a;eu\U{oo,t)\n), 

where the path of the integration in the complex E plane is given in Fig. 4. 
The energy dependent function V2a,n{,E) is the vertex function for the decay 
of the shell model states | n) into the two a particle states and the energy 
dependent function Rn{E) is the self energy of the state \n). They are all 

^■"'^It is possible that some states in the continuum also have large matrix elements. In 
case of the A=8 system, it does not turn out to be the case. It will be demonstrated in 
the following that the direct transition from ^B or ^Li to two a particle scattering states 
characterized by a parameter C is small. 

^■^i.e. the energies at which the phase shifts of the a-a scattering S-matrix pass through 
{n + ^)tt from below, with n = 0, ±1, . . .. 



53 



related to AV in the following way 

R^{E) = {n\AVF\n), (D.6) 

T2a,n{E) = {2a\AVF\n), (D.7) 

where the state F\n) =\ ipn) satisfies a Lippmann-Schwinger type of equation 

l^n) = K)+ p rr^ AV]^^), (D.8) 

with I n) an eigenstate of H shell- 
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Figure Captions 



Figure 1: Diagrammatic representation of tlie /9-decay in A=8 system. 

Figure 2: ^Be energy levels calculated in the shell model compared with the 
experimental ones. The shell model energy level at 13.37 MeV should be 
identified with the experimental one at 16.92 MeV. It is explained in the 
text. The shell model energy level at 16.19 MeV should be identified with 
the experimental one at 16.63 MeV. 

Figure 3: Experimental resonant levels in the ^Be* system with J^ = 0"^, 2"*" 
and 4"*". The first column is the O"*" series, the second column is the 2~^ series 
and the third column is the 4"*" series. 

Figure 4: Result of the least square fit of the R-matrix parameters to the 
experimental phase shifts of the J^ = 0^ state. 

Figure 5: Result of the least square fit of the R-matrix parameters to the 
experimental phase shifts of the J^ = 2+ state. 

Figure 6: Result of the least square fit of the R-matrix parameters to the 
experimental phase shifts of the J^ = 4~^ state. 

Figure 7: R-matrix calculation of the a particle spectrum in the /9-decay 
process ^B -^ a + a + e"*" + z/g- Solid curve represents the full calculation and 
dashed curve represents the calculation without the final potential scattering 
between the final two a particles. Open circles represent the experimental 
measured values. The contributions of the 3 MeV state and the 16 MeV 
doublet states to the total counts are also indicated. 

Figure 8: R-matrix calculation of the a particle spectrum in the /9-decay 
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process ^Li -^ a + a + e + Ve- Others are the same as Fig. 8. 

Figure 9: R-matrix calculation of the a particle spectrum in the /9-decay 
process ^ B -^ a + a + e'^ + u^. A linear scale is used. Open circles represent 
the experimental measured values. 

Figure 10: R-matrix calculation of the a particle spectrum in the /3-decay 
process ^Li — > a + a + e~ + z/g- Others are the same as Fig. 10. 

Figure 11: Energy dependence of ?76(2''', 0+). 

Figure 12: Energy dependence of 779(2+, 0+). 

Figure 13: Energy dependence of r76(2+,4+). 

Figure 14: Energy dependence of 779(2+, 4+). 

Figure 15: 3-dimensional plot of Ri'^^^^{Er,ee) for the /3+ decay of A=8 
system. The vertical scale is 10^^. 

Figure 16: 3-dimensional plot of R2'^^^^{Er,ee) for the (3^ decay of A=8 
system. The vertical scale is 10^^ 

Figure 17: 3-dimensional plot of i?f^^^^(i?j., Ce) for the f3~ decay of A=8 
system. The vertical scale is 10~^. 

Figure 18: 3-dimensional plot of R2^^^^{Er,ee) for the l3~ decay of A=8 
system. The vertical scale is 10~^. 

Figure 19: Averaged PTRV response function R^'^^^ (Er) as an function of 
Er. Solid lines represent Rg + Rn . Dashed lines represent ^9 and Rn for 
the P~ decay. Dash-dotted hne represent Rg and R\i for the /S"*" decay. 
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Figure 20: Averaged PTRV response function R2'^^^ [Er) as an function of 
Er- Solid lines represent Riq + R12 . Dashed lines represent R[q and R12 for 
the (3~ decay. Dash-dotted line represent R\q and R12 for the (3~^ decay. 

Figure 21: 3-dimensional plot of Rf^i^{Er,ee) for the /S"*" decay of A=8 
system. The vertical scale is 10"'^. 

Figure 22: 3-dimensional plot of -R^f^(i?r, Ce) for the P~ decay of A=8 
system. The vertical scale is 10"*^. 

Figure 23: 3-dimensional plot of i?^f^(_Er, ^e) for the P~^ decay of A=8 
system. The vertical scale is 10^^. 

Figure 24: 3-dimensional plot of R^^2^{Er,ee) for the P~ decay of A=8 
system. The vertical scale is 10^'^. 

Figure 25: The path for the complex E integration. Solid dots represent poles 
of the integrand and the dashed line represents the cut of the integrand. (It 
is usually located on the real axis extending all the way to positive infinity. 
Two cuts along different directions are equivalent provided that there are no 
singularities in between them.) 

Figure 26: Diagrammatic representation of the resonant level contribution 
to the a-a scattering amplitude. 
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Table 1: Kinematic functions considered. Here "T" and "P" denote the time 
reversal and parity reflection transformations. 



fT = 


^\r 
p^-^ 


Even 


Odd 


Even 


K,,{0) K^t\01) KUOOO) KUOll) 


X 


Odd 


X 


X 


cr = 2 


^\r 
^^^ 


Even 


Odd 


Even 


KU2) K^Pi21) K^sv\22) X„„(211) 
X™(202) i^™(220) X™(222) 


X 


Odd 


X 


xi+)(22) Ki-\2l) X™(212) i^™(221) 
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Table 2: Expansion of the reduced matrix elements of the non-relativistic 
hadronic one body weak current multipole operators in terms of power series 
of the momentum transfer k, = |ke + k,^|. Here ^^ = fi'A ± oXf fi'^- 



Vector Current 



Axial Vector Current 



Cn 



^^^^-Ijp^\ 



CI 



K 



'^^\-^^^9AT9T)A, + \g^A, 



U 






L\ 



T^'' y"' 



6M^ 



rAv 






A, 



Ml 



IK 

2M' 



G 



M] 



'2 , ifi ^ . 

-A2 F1A4 

3 3M 



E! 




QM' 



rAv 



V2k^ 



A, 



Co 



1 k' 
T5M- 



:FiA, 



Ml 



1 k' . 



1 kA 
15 iW^ 



FiAe 






2i k' 

zdA^rpjAv) 



15v/7""M^ 



Eo 



1 /3 kA 
T5V2M^^^^^ 



-^3 



T5V7^^M^^^° 
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Table 3: Definitions for tlie static multipole operators. Here (J/ ||. . .|| Jj) de- 
notes tlie reduced matrix elements of the corresponding multipole operators. 
Vector 0" represents the collection of three Pauli 2x2 matrices. 



Al 


y,{Jf\\r^^^Y,\\J,) 

i=l 


Ae 




A2 




A, 




A, 


y^Jf rWrYla. VII J.) 
1=1 


As 


M^y{Jf\\r^^\^y,^.a J,) 


^4 


y,{Jf |r(±V3^n-V| J.) 


A, 


M2y^(J^||rWr23^22-a J,) 


A, 


M'y{Jf\\r^^h-'Yo J.) 


^10 
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Table 4: R-matrix parameters for the J^ = O'^ state obtained from a best 
fit to the experimental a-a scattering phase shifts. 



a MeV-'^ b MeV-^ c MeV-^ wq (MeV) tq (fm) z MeV-^ 
n QO/1 n 070 K «oo n ni i /i Qn i Q/1 x^ i n-4 



E (MeV) 



0.0 



0.324 



-0.272 



-6.622 



0.011 



4.39 



1.34x10 
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Table 5: R-matrix parameters for the J^ = 2+ state obtained from a best fit 
to the experimental a-a scattering phase shifts. No energy dependence for 
the doublet states located near _E ~ 16 MeV is considered, which is indicated 
by O's in the table. The a-a phase shifts are very insensitive to the values of 
Wq for the 16 MeV doublet states in the energy region considered, they will 
be determined in the a particle spectra fit. 



a MeV-^ b MeV^^ c MeV'^ wq (MeV) ro (fm) z MeV'^ 



E MeV 



3.36 


0.145 


-0.01704 


0.00216 


1.80 


16.63 











? 


16.92 











? 



2.45 5.2x10"^ 
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Table 6: R-matrix parameters for the J^ = 4"^ state obtained from a best 
fit to the experimental a-a scattering phase shifts. 



a MeV~'^ b MeV'^ c MeV'^ wq (MeV) ro (fm) z MeV- 



E MeV 
12.15 



-0.400 



-0.048 



-0.126 



3.54 



3.65 



0.0 
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Table 7: Values for various static single nucleon weak current form factors 
used. 



-1.254 1.00 4.71 0.00 
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Table 8: Kinematic functions correspond to the contribution of the symmetric 
(with respect to A^^ and B^^") and dispersive parts of the final lepton-hadron 
Coulomb scattering T-matrix contribution. 



cr = 2 




Odd 


Even 


Odd 


Kssi2) Ki+\21) Ki-\22) K^4211) 
if™ (202) i^,,(220) if™ (222) 


X 


Even 


X 


Kit\22) k!,-\21) X„„(212) if,„(221) 



